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Modeling of General Surface Junctions of Composite Objects 
IN AN SIE/MoM Formulation 

Joon Shin, Allen W. Glisson, and Ahmed A. Kishk 
Department of Electrical Engineering 
University of Mississippi, University, MS 38677 

Arstract This paper discusses the modeling of various kinds of surface junctions in an 
SIE/MoM (Surface Integral Equation / Method of Moments) formulation appHed to complex 
objects consisting of arbitrarily shaped conducting and dielectric bodies. Methods of describ¬ 
ing various types of junctions and systematically incorporating them in numerical solutions are 
presented. The procedures are of interest for the specific application of arbitrarily shaped di¬ 
electric resonator antennas and their associated feed structures and packaging. An E-PMCHW 
formulation in conjunction with a moment method procedure using generalized triangular basis 
functions.^ presented to deal with such general junctions. 

I. Introduction 

The modeling of general surface junctions in an SIE/MoM formulation is considered in this 
work. The specific application leading to this study is that of a dielectric resonator (DR) antenna. 
Since an experimental study of a cylindrical aelectric resonator (DR) antenna was reported in 
1983 [1], this antenna has drawn continued interest because of its small size, efficiency, and 
potential ability to perform multiple antenna tasks via simple mode coupling mechanisms, ^e 
configuration of a DR antenna may range from a very simple one which allows analytic solutions 
to a very complex one. A typical structure for a DR antenna is a DR element of high ffielectric 
constant excited by a single feed such as a microstripline or coaxial cable. Various shapes and 
combinations of DR elements as well as various feed structures have been suggested, however, 
which may improve the antenna performance in the areas of bandwidth, power handling, and 
antenna efficiency. 

Rigorous SIE analysis methods for non-trivial DR antenna configurations have been available 
mainly for body of revolution (BOR) objects [2,3]. This work results ffom an int^est in the 
analysis of DR antennas of more arbitrary configurations, which may include genei^ 3D objects 
comprising an arbitrary combination of conducting and/or dielectric bodi^ of arbitrary shapes, 
nsiTig an SIE/MoM method with triangular patch basis functions. The junction modeling problem 
has been considered in previous works for conducting surfaces [4], for ffielectric surfaces [5], simple 
combinations of both for BOR objects [2,3,6,7], and for a general case of conducting, diele^nc, 
resistive, and impedance boundary conffition surfaces [8]. This work attempts to provide a 
formalism for systematically describing junction models for a wide variety of junction lypes. 


II. Formulation 

A. Problem Description 

The geometry under consideration is a general inhomogeneous body with Nr dielectric re^ons, 
each of which may contain conducting bodies as weU as impressed sources as shown m Fig. 1. 
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The regionshavepennittivitiesf..andpera,eabmties//.., where.-= 1 ar 

may be complex to represent lossv materifl.]<! Mnn i i ^ ^ 
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y # 0; however, the normal unit vectors n, and fi, are in opposite^^rior^ e^ 
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Fig. 1. General geometry under consideration. 
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B. The Field Equivalences 

According to the equivalence principle [9], the original problem can be decomposed into Nr 
auxiliary problems, one for each dielectric region. To obtain the auxiliary problem for region Hi, 
the impressed sources of the original problem are retained only in region Hi and the bound^es 
of the region are replaced by equivalent surface currents radiating in a homogeneous medium 
with the constitutive parameters of region Ri. Electric curr^ts are used for the conducting 
surfaces while electric and magnetic currents are used for the dielectric boundaries. The electnc 
and mametic currents appearing on opposite sides of a dielectric interface in differrat auxihaj^ 
problems are taken equal in magnitude and opposite in direction to assure the contmmty of the 
tangential components on these boundaries as they are continuous in the onginal problem. In 
this procedure, the fields produced within the region boundaries by the equivalent currmts and 
the impressed sources in region Ri must be the same as those in the original problem, whde the 
zero field is produced outside these boundaries. The electric and magnetic currents along 5,- are 
m X Hi and Mi = Ei X rii, respectively. . . u j 

A system of surface integro-differential equations can be obtained by enforcmg the boundary 
conditions of continuity of the tangential components of electric field on the cpnductmg 
and both electric and magnetic fields on the dielectric surfaces. This results m the &PMCHW 
formulation [61 when there is no junction in the problem. For problems havmg g^eral junctions, 
however it is not easy to express the integral equation system exphcitly apa^ from tte t^tmg 
procedure. Thus the system of integral equations is presented in the next section after descnbmg 
the junction modeling and the basis functions. 


a Modeling of Junctions in the Moment Method Solution 
Arbitrarily shaped surfaces are discretized in triangular patches and the equivalent surface 
currents are approximated by expansions in the RWG basis functions on the patches, which are 
expressed as [10] ^ 

( 1 ) 

T»=l 






T G jS^ 
otherwise, 


is the Diunber of electric basis functions, and are the positive/negative dom^ or the 
rem-Zto- faces of the basis function, respectivdy. For magnetic currents, I* 

ieBned similarly. The testing functions tV and TJ- are also taken to be the same as (2). With 
the basis and testing functions defined we have a matrix equation 


TyrrH r I/O 

'TmT„ j/Tsn 


When there are general surface junctions, the current related to an unknown coeffiaert may 
idst on many diffment surfaces. In such cases, the expression (1) “ ^ 

jcample, there is an electric current on a adectric surface m the repon ft eqmval^ 
ad mother one flowing in the opposite direction m the region ft problem, represented by » 
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as shown in Fig. 2(a). The expression in (1) for the electric currents has this sort of implication 
for the basis functions when the dommn of the unknown involves a dielectric interface. When 
more than two dielectric surfaces meet at a junction, this scheme does not work. Thus for general 
junctions, we seek another way of expressing the generalized current more rigorously. We will 
use two different basis fonctions for the same unknown coefficient related to a dielectric surface 
as shown in Fig. 2(b). In other words, the unknown co^dent has a multiplidty of two when it 
represents the electric or magnetic current on the dielectric face. Extending this to the general 
case, the generalized current is defined in terms of the generalized basis functions as 

C(r) = {J(r),M(r)} = {§/„B„(r), f) /„B„(r)} (4) 

n=l n=l+NTj 

where, 

B„(r) = B*^(r), B,^(r) = B2(r), Tn = r*, with ft = n, 

Bn(r) = Bj’”(r), B,^(r) = Bg‘(r), r„ = r*, withft = n-iVrp 


if n<NT^ (5) 

if n > Ntj (6) 


N = Ntj+Nt„ 
B?W = 

V=1 


r=l 


B2.B& 


rhf 


Rn, 


= the V** basis function of u = 1,..., Tn 

= RWG basis function defined over the corresponding patches as in (2) 

= multiplicity of the unknown coefficient, J„ — 

= total number of faces (surfaces) connected to the junction for /„ 

= number of dielectric faces (surfaces) related to In 
— firom-face and to-face of B„^. 

= region of B^. 


= “’hf 

n^n +otherwise 


(7) 

( 8 ) 








(a) Conventional representation 


(b) New representation 


Fig. 2. Two methods of representing ba^ fonctions. 
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Notice that there is one-to-one correspondence between or and the parameter set 
{//n,>i/nvj-Rrh,}- Th® niimhprfi of unknowBS and basis functions for a given junction or edge 
are deterained from the types and numbers of the faces connected to the junction by considering 
proper boundary conditions at the junction. The methods of determining them and systematical¬ 
ly incorporating them in the MoM solutions have been developed. Examples of modeling general 
junctions are shown in Fig. 3, where Jn and Mn are used instead of Bn^ and B^"*, respective- 

ly. Tte generalized testing functions {t5}Si. and are ^ defined in a 

similar manner. We also define Q, the generalized current for the region Ri equivalent problem, 

C((r) = {Ji(r).Mi(r)} (9) 

where 

JiW = E (10) 

n=:l »=1 

Mi(r) = E /„E«tiB^(r;//»..t/».,lJ».) (11) 

n=://7^.+l V^l 

{ 1 Rfiv “ Ri 

o’ Sherwise 

With the set of basis functions in (4-8), one may apply the boundary conditions of tangential 
field continuity at each sub-domain of the basis functions. By merely applying the boundary 
conditions, however, the total number of equations may be greater than the number of the 
unknowns because of the multiplicity of some unknowns related to junctions. The usual methods 
of solving equations apply only when the number of equations equals to the number of unknowns, 
N. Such a set of N equations can be obtained by taking the n*^ integral equation as the set 
of simultaneous integral equations (or summation of them) which satisfy the proper boundary 
conditions on the subdommns of the basis functions (Bm,,u = l,...,r„) related to the unknoira 
coefficient, J„. It is possible to get such a surface integral equation system by testing with the 
generalized testing functions as follows 


No Tm ' « • « 

5](E?(Ci)|tan>ETm,, ) = -2(^!tan» >1 = 


5:(n,xi^(Q). = 


i=l «=1 

m = Ntj +1, Ntj + Nt^ , 


( 12 ) 

(13) 


where, 

(tg) = 

f 1 Ri = Rr^ 

S^i = field contribution coefficient = | otherwise 

= region of the testing function, Tm,, 
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Fig. 3. Modeling of general surfece junctions. (Ct,» = 1,2,3,..., is an entry-counting index.) 






Fig. 4. A test case, (a) PEC alone (b) with phantom dielectric. 


and Hf) and (e|, Hj) are the scattered fields due to C,- and incident fields, respective¬ 
ly EquaW (12) and (13) are the E-PMCHW formulation [6] extended to general junctions. 
Substituting C,- into (12) and (13), the impedance matrix and excitation vector elements m (3), 
ZS? and vS', for example, are expressed as 

^ ([Ef(^ (r^ //n„,tfn,,i2nv))!tan> 

n = l,...,iVT^- and m= l,...,iVT^ (14) 

i=l V=1 

Some subroutines of EMPACK [11] have been used for the integrations over the triangular do¬ 
mains which appear in (14) implicitly. 


III. Numerical Results 

A T-shape junction of three 0.1-m wide and 0.3-m long PEC strips is tah^ as an example. For 
comparison a semi-circular cylinder of phantom dielectric having 0.1-m hei^t and 0.3-m radius 
is attached to the T-shape junction as shown in Fig. 4. The z-directed surface currents along 
the contour lines, (-0.3,0,0)-> (0.3,0,0) and (0,0,0)(0, -0.3,0), located at theater o^ 
strip are computed for a plane wave excitation. The plane wave expressed ^ ’ 

where jt* = -^cos(^*sin^-j/sin^’sin0<-zcos^, Eo=F?(icos^cos-;6*+ycos^sm<^-zsm^), 
^ = ^< = 45% £5 = 1, A:p = 27r/v^, and / = 300 MHz. The results in Fig. 5 show very good 
agreement as well as the expected singularities at the end of the strips. 

IV. Conclusion 

A systematic procedure for modeling of the general junctions of any combination of condu^g 
and/or dielectric bodies in an SIE/MoM formulation has been presented. With the succe^ 
modeling of general junctions, it is possible to ^ply the E-PMCHW formulation to a large class 
of problems induding dielectric resonator antennas of complex configuration. 
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Fig. 5. jf-^ected current denies along the contours (o ; PEC alone, x ; with phantom dielectric). The 
arrows denote the start of the second contour. > h 
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